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Technical Note 



Vessel Enhancement Filtering in Three- 
Dimensional MR Angiograms Using Long- 
Range Signal Correlation 



YipJng P. Du, PhD • Dennis L. Parker, PhD 



Sttl*n veiwcia In three ^Iffl^n- 
■Lonal MR- angiograms have 
low TUibtUt7 ta n*«rtfiM»n4n- 
tcniltj projection Images fce- 
cvuse of their Iitw contrast, h* 
a previous study, we had two 
nonlinear Wets that ap- 
peare4 to give significant im- 
provement In small vessel 
detail. In this paper, we report 
on a generalization of this fll^ 
t*r that allows a mote general 
modeling of the vosacls and a 
mox* complete suppression of 
bac&gjfnund. One implcmefl" 
tatioA of the general filter 
gave a Teasel mean eontraat- 
to-Aolse ratio that t» 2.tf2 and 
3.51 time* higher than the 
veuel mean cofltraBt-to^noise 
ratio obtained nalng our pre- 
vionaty reported maximnm* 
minimum (nmx-mln) filter and 
cross- taction filter, respec- 
tively. 

|#i£*taznic Cocitre* fmhanocmenc • ttn-* 
age Eftertng * tma^ (tfixeatsUi^ • V»flCulwr 
MJ&t* ♦ InterpbtetlOC 

JHBJ 1SOT; 7:447-*SO 

AUtwertiUc**: MRA. = MR sn£qcmm. 
TOF - tttne-cflUght MP » m«*imim»-BV 
tuwtfy prrjectkm. CHft » cotrtrw^to- 
nattt ratio, 2PI - Mfo-fifled interpolation. 
5D - second dl3cT«U* 
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THIS PAPER IS A CONTINUATION of our 
studies to use Image processing tech- 
niques to improve small vessel visibility tn 
magnetic resonance angiography (MRA). 
In tfmc-oMVght fTOF) MBA techniques, 
vessel signal Is low In small vessels, where 
the slowly flowing blood experiences many 
excitation pulses* The vessel- background 
contrast of small vessels is fujther re- 
duced by the use of the maxinium-rnten- 
s«y projection (MIP) algorithm (1-4). 

Limited improvements In vessel visibility 
have been, obtained by using different dig- 
ital filter* (5-7) and other postprocessing 
techniques (B-10). The development of 
more effective alters will be necessary be- 
fore toe filtering technique can be used to 
improve clinical diagnoses, 

We have recently reported on two non- 
linear, anisotropic, digital vessel-en- 
hancement filters (11,12) that were 
developed based upon the characteristics 
of small vessels In a nearly uniform back' 
ground Signal distribution and upon the 
properties of a second difference opera- 
tor. The filters were designed based on a 
line- segment model in which a small ves- 
sel Is approximated as a sequence of vox- 
do tfe. a line Geg)cne*tt) with higher signal 
Intensity than the surrounding back- 
ground. Although these filters were 
shown to improve significantly the visi- 
bility of small vessels, both filters tend to 
increase image noise. In Bpite of the ef- 
forts that were made to suppress noise in 
the design of the niters. 

tn this paper, we employ a more general 
vessel model that approximates the vessel 
as a cylinder segment of finite width. This 
model can be used to design a family of 
filters that are optimized to enhance ves- 
sels of sped fie sizes and suppress noise 
and signal from stationary tissue. We 
demonstrate that one of our previous fil- 
ters is a nonopttmized specific example 
from this family of flltere. Based upon this 
generalized model and our knowledge of 
the image potDt spread function, we have 
developed a more optimal nonlinear, ani- 
sotropic vessel-enhancement Biter that 
uses bog-range signal correlation and is 
designed to enhance vessels of the small- 
est diameters thai, can be observed tn our 
MRA images. Quantitative contrast-to- 
noise-ratta (CNft) measure avtnU have 
been used to assess the improvement of 
vusel visibility obtained by using the pro- 



posed filter compared with our previously 
reported HlrerB. 

• THEORY 

Modeling "j^uve-DbneneUmal AWA 

Out- generalized filter Is based upon the 
signal properties in the tfcuree-dimenaianal 
TOF MRA images, where flowing blood ex- 
periences fewer RF excitation pulses and 
therefore has higher signal Intensity than 
the *u!T0Un4ing background that IS highly 
saturated and exhibits very little spatial 
variation. Image noise consists of ran- 
domly distributed Intensity variations over 
the entire image- Furthermore, the filter 
also takes advantage of the improved CNR. 
continuity and smoothness of small ves- 
sels that results from zcro*fllled interpo- 
lation [ZFll (13-16). We believe that ZFI of 
MRAs is a critical step to ensure the suc- 
cess of subsequent vessel enhancement 
Altering techniques (11.12). 

In ZFt interpolated image data sets, me 
image point spread Function is about the 
same width as the interpolation factor, 
and the crOss^section of a small vessel is 
represented by multtple grid points. There- 
fare, the line model ta Figure 1 does not 
adequately represent a ornafl vessel when 
the thr**-dImenelon&l angiograms have 
been interpolated before the filtering pro- 
cess (12). We model a small vessel as «. cyl- 
inder that contains voxels with higher 
intensity than barAgroiirrd. In a plane con- 
taining the vessel cross-section, the voxels 
in the center of the cylinder havetne high- 
est Signal Intensity compared with the vox- 
els dose to the edge of (but anil within) the 
cylinder. We model the stationary tissue 
as a constant background and assume 
that the correlation of tmage noise is neg- 
ligible from voxel to voxel. 

In our previous ZFI studies <i?,l£). we 
used the Interpolation factor. L. to de- 
scribe the power of magnification in the 
matrix dimension. For instance, the inter- 
polation from a 236 X 236 X 32 to 1024 
X 1024 X 128 has L <= 4- We found that 
in a three-dimensional MRA interpolated 
using ZFI With an Interpolation factor of U 
the apparent diameter of the smallest vis- 
ible vessels is about 2L-1 pixels. There- 
fore, this cylinder model would appear to 
be correct jar even the smallest vessels in 
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Figure 1. Schematic shows the orienta- 
tion correlation of a small vessel between 
three neighboring 3x3x3 Cubes (fC = 
1). Each of the arrows represents the di- 
rection tn which hap the maximum 
value in this 3x3x3 cube. 

the interpolated MR angiograms. Qrven 
this observation, we would propose that 
an efficient filter 1 (One that gives good Ves- 
sel enhancement and soft tissue suppres- 
sion) win approximate the properties of the 
desired vessel limited by the imaging point 
spread function. 

A Simple Generalized Filter 

To test our model* we have implemented 
an approximation of the desired Biter for 
which the coefficients conajst of three pos- 
itive central points and two negative points 
at an appropriate distance from the cen- 
tral points. This vessel enhancement filter, 
which is applied to the three-dimensional 
image data before MIP reconstruction, 
uUUzcs the 13 discrete directions that 
pass through the center of a three -dimen- 
sional cube with 3 voxels on each side (sec 
second figure in Du et aJ. (12]). These 13 
discrete directions can be written in Car- 
tesian coordinate as: [1, 0, 0), (0. 1. 0)„ (0. 
O. 1), (t. X. 0). (1. - J, O). (1. 0. 1). (1. 0. 
-I), (0. I. 1). (0. -1. 1). (1. 1. X). (1. -t. 
1). (1. 1. -1), and (I. -1. -X). The algo- 
rithm begins by applying a long-range dif- 
ferential operator, D M , along each of these 
13 discrete directions. The D w operator at 
an image point t and operating along a dis- 
crete direction J la defined as. 

where f f ts the image Intensity at point t, f, 
_ , and J,+ , arc the image intensities at the 
neighboring points of J, along direction 
j. and Mis the signal correlation range of 
the operator. The second difference (SD1 
operator, used In Du et at ( 12) Is a special 
case of D M when Af = I . 

Based on the cylinder model described 
above, the operator haa the following 
properties, which arc essentially identical 

to those considered in our previous filter 
design 02); (a) Ai a point centered on a 
vessel, the O m operator has a large poslLtve 
value tn a direction perpendicular to the 
vessel (j„ in Fig. 2J and a small positive or 
negative value m ihe direction parallel to 
the vessel \J b in Fig. 21; 2) At a point near 



the vessel edge, the O u operator has a 
large negative value tn a direction perpen- 
dicular to the vessel edge (y c in Fig. 2) and 
a small positive or negative value tn the dl» 
rection parallel to the vessel (j a in Fig. 2); 
3) Tn background regions with small inten- 
sity Variations, the D„ operator has a small 
value, either positive or negative, in aXi di- 
rections: 4) Finally, the r\, operator has 
large positive values in all directions on 
positive noise spikes and large negative 
values tn all directions on negative noise 
spikes. Using these properties of the 
operator, we have designed a family of fil- 
ters (with different correlation ranges) to 
reduce image background, enhance vessel 
structures, and suppress noise spikes. 

The- Oy operator enhances cylindrical 
structures with diameter of f2M-l) grid 
points or smaller. The D u operator is then 
used to construct a g M operator, similar to 
using the SD operator to construct a 
cross-seotlon filter: 



1 * 

9vm = - 2* d 8uk ~ mew 10.1)^1 



(2) 



where N is the number of discrete direc- 
tions that are perpendicular to j. and 
Oyji is denoted to the operator at point 
i along a direction perpendicular to j. 

the oaximum value of g uW over the 13 
discrete directions provides an estimate of ' 
whether the image voxel ts vessel or nofi- . 
vessel. If this image voxel Is vessel, the 
vessel is most likely In the direction j at 
which has the maximum value (sec 
Fig. W Instead of using to construct a 
filter similar to the Cross^ectlon filter, we 
designed a filter using the relationship be- 
tween g lhM and the orientation of vessels. 
In the design of the filter, we also assumed 
that the small vessel has a strong OTlen- 
taUon-dependent correlation within a 
small cubic volume of [2{M + K) -1] x 
U(M+K) -i| x t2(M-ffC)-l| voxels, The 
Alter with signal correlation range of Wand 
orientation-dependent correlation range of 
K is given by: 



„ - ~ K * 



where L "(i + M is the summation 

of Pyj, along the direction of J. 

For the visualization of the smallest ves- 
sels In a thxee-dtmens lo nal MRA Image 
data set with an interpolation factor of L. 
we tested the values of M => L and K = I. 
- 1. These values arise because or the ex- 
pected correlation distance of 1L — \ that 
arises from the interpolation process. The 
chosen values allow averaging over the ex- 
pected correlation distance and the SD to 
be taken over a distance that ts just 
greater. 

* METHODS 

To test the vessel enhancement filtering 
techniques, a three-dimensional single 
slab TOF MBA image was acquired on a 
conventional GE Signa MR Scanner (CE- 
Mcdlcal Systems. Milwaukee. Wl) that op*, 
crates at 1.3 Tesla. To reduce echo time, 




Figure z. schematic shows the location 
and direction that the differential operator 
Is applied, j j on the vessel centcrlinc 
perpendicular to the vessel; on the ves- 
sel centerltne, along the vessel: on the 
edge of M the vessel, perpendicular to the 
vessel; J*: on the vessel centerltne. along 
the vessel. 



the readout was made asymmetric auch 
that the echo center was placed at a point 
before the center of the readout window. 
For the studies of this paper.- the shift was 
25% of the total readout time. The circle of 
WUUs was imaged In a 33-yesx-old human 
volunteer with the three-dimensional TOF . h , 
technique (25B X 256 x 32 acquisition ^ '■ ' f 
matrix. TfUnTsed/TEfrnsec) - 46/7.0, 25* : ■'. 
flip angle. 18-cin ftcld of view, 0.75-mm 
Section thickness). 

In this study, we applied ZFI to increase 
the Image matrix dimensions in all three 
directions by s factor of two fee. & - 2) fox 
the image data before filtering. After ZFl. 
three sets of filtered images were obtained 
by using the previously reported maxi- 
mum-minimum and cross-section Alters 
(11,12) and the new Cfi flutter with Af = 2 
and K = l. The resulting image data were 
displayed using the MIP algorithm. 

A vessel CNR measurement technique 
(1°) was used to assess the CNR at each 
point along a vessel segment and the mean 
CNR of this vessel segment tn these three 
dimensional filtered data sets. Vessel con- 
trast is the difference tn image intensity 
between a vessel point and its surround- 
ing norrvesscl background.. The standard 
deviation of image intensity In a three-di- 
mensional region of background was used 
as the estimation of image noise for these 
three-dimensional filtered Image data sets. 

The Cy filter has been Implemented in 
the C language on a Sun" Sparc 20 work- 
station (Sun Microsystem*. Mountain 
View. CA). The computation time was 20 
minutes for a 5 12 X 512 x 64 data set. 



• RESULTS 

The Images of Figure 3 are al) obtained 
from the same three-dimensional MRA im- 
age data. Figure 3a shows the MIP of the 
original 512 x Sl2 x 64 interpolated Im- 
age. Figures 3b. c and d show the MIP Im- 
ages of the three-dimensional filtered 
image data sets using ihe marimunVmin- 
Imum filler- cross-section filter, and the 
Ci filter, respectively. Although Inert is 
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Hgure 3. (a) The axial MIP image of a fcero-ffllcd -interpolated thiee^dimenslonal MRA 
Image data act: (b) The MIP Image of the filtered three-dimensional MftA data using a 
maadmunHttitntmum fliter (e) The MJDP Image of the filtered three-dimensional MBA 
data using a craaa-acctum, filter; (d) The Mip image of the filtered thrcc-cUmjensional 
MRA data using the G J filter. The same threc-dWnelonal MRA data as In W was used 
as the original linage data for these three fijtcna. 



" v some variability In small vessel visibility 
between the imagea, It is evident thAt, In 
general, the small vessel detail ts better 
and more continuous in Ftgure 3d. Figure 
3 t* a typical comparison of the effective- 
ness of these three Biters In suppressing 
background signal and noise. In Figure 3. 
wc found that the G± Alter considerably 
Improved the srncothjiesa of vessels and 
reduced the noise level. We also found that 
the signal toss hi large Vessels was rtotice- 
ably reduced by using the filter 

com- 
pared with that obtained using the fnaxH- 
mum-rofcrtfTnum and cross- section filters. 

Figure 4 plots the image intensity profile 
along a horizontal straight line across the 
posterior ccniro unseating arturJea in the 
MIP images (as shown in Fig, 3a). Hie 
maximum Intensity In each of the profile 
plots was normalized to 100. The Intensity 
profiles on JMJP Images of the Altered three- 
dimensional image data; seta show the ves- 
sel enhancement by using the maximum- 
minimum filter and the GA filter. Small 
vessels labeled as G, H, and become vis- 
ible In the MIP linage of the maxlmum- 
romlnnim filtered data, to Ftgure 4a, we 
- observe that the background signal waa 
well suppressed. However, the noise level 
in background tissue xuid fca air wto in- 
creased. The plots to Figure 4b show that 
the After suppresses background signal 
and noise more effectively compared wtth 
the maxl mum-mtaiinufii filter. The vessels 
labeled as C. G, H. and K have Increased 
vessel-background contrast In the Image 
processed using the CJ filter compared 
with the image processed using the maxi- 
mum- minimum filter. 

Figure 5 Illustrates the comparison of 
CNR meaaurementa for three filtered im- 
age data sets. All of the vessels for which 
CNR could be conveniently measured were 
divided into 53 -segments. These segments 
included large vessels, such as middle ce- 
rebral artertea, and the vessels barely visible 
la the local maxbmmHntensity projections. 
The mean CNR measurement obtained 
from these S3 vessel segments In the max- 
imum-minimum filtered Image data was 
plotted as the xaxis. The mean CNR meas- 
urements from the same vessel segments 
in the image data sets processed using the 
cross-section filter and the CJ filter *rerc 
plotted along the y sods, respectively* In 
this figure, large vessel segments usually 
have high CNR values, and small vessel 
segments usually have low CNR values. 
The Ieaat-meanraquare fit obtained using 
a linear function, y = ox was performed 
for the mean CNR measurements, it was 
found from Figure 5 that the MRA pro- 
cessed using the Cj Alter' has a consider- 
able improvement in CNR compared with 
the same MR angiogram processed using 
the previously reported filters. In this 
study* the vessel 'mean CNR obtained us- 
ing the Ci filter was 2.53 and 3-51 ttmes 
higher than the vessel mean CNR obtained 
using the me^dmunv minimum filter and 
cross -section filter, respectively. 

• DISCUSSION 

The results presented here substantiate 
our prediction that small vessel visibility 
will be improved by using a filler thai 
matches the expected signal of the small 



vessels In the Interpolated Image data set 
Although, as observed above, some small 
vessels may appear more visible In other 
processed images, in general, the simple 
implementation of our generalized fitter 
(the C£ filter) waa found to considerabfy in- 
crease the visibility oJt small vessels and to 
suppress stationary background signal 
and Image rtol&e more effectively compared 
wtth our previously reported filters (12). 

Our modeling of small vessels in three* 
dimensional MRA data is based on the flow 
pattern of blood. The blood flow In a small 
vessel is usually fairly stable, so that this 
flow can be assumed to be laminar. In 
laminar flow, the fluid velocity has a par- 
abolic distribution, wtth the maximum ve- 
locity at the center of the vessel and zero 
velocity at the edge of the vessel. In three 
dimensional TOF MRA imaging tech- 
niques, faster flow experiences fewer ra- 
dJofrequency excitations and has less 
signal saturation. Therefore, the signal 
contribution from the central part of the 
Row is higher than from the edge part of 
the flow, fa the original tmage data sets, 
the image grid points are too coarse to 
show Such Signal distribution In small 
vessels. With 2F1, small vessels have a 
considerably Improved CNR, continuity 
and smoothness. In a three-dimensional 
MRA wilh ZFt, the cras^-scctton of a small 
vessel is represented by multiple voxels. 
After ZFI. especially when the interpola- 
tion factor is large, the signal distribution 



of a small vessel is much smoother tn all 
directions compared with vessel smooth- 
ness In the original Images. Realistic mod- 
eling of a small vessel in such Interpolated 
images should treat a small vessel as a cyl- 
inder instead of a Unc that conefots a se- 
quence of high Signal voxels. This cylinder 
has higher signal intensity along the cen- 
ter tine and lower signal intensity close to 
the edge. For a small vessel, the diameter 
of this cylinder in number of grid points is 
approximately the same as the Interpola- 
tion factor. The long range D u operator was 
derived from this modeling of a small ves~ 
sel. 

As an example, a hypothetical signal in- 
tensity profile of a small vessel across the 
cross-Section of the vessel is plotted In Fig- 
ure 6a. At the center of this vesael. the SD 
operator has a value of 14 [— 2 X 130 - 
123 - 123). The fJj opera tor defined in Eq- 
(11 has a value of 60 (- % x 130 - 100 - 
1001 at the same location. This example 
shows the improved capability of enhanc- 
ing small vessels using the G£ filter. 

In our prevtous studies, we observed 
that large vessel signal waa decreased in 
general when the maximum-minimum 
and cross-section filters were used. In the 
current study, twe also observed signal loss 
in large vessels, but this loss was appar- 
ently less than that from the other filters. 
In the eg filter, the Du operator is sensitive 
to cylindrical structures with diameter of 
2M - 1- grid points, whereas the 5D oper- 
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Figure 4. The comparison plots of troagc 
Intensity profiled for the M1P. Images of (■) 
original Image data and the maxbnum- 
mJnimwm Altered Image data: (b) the max- 
famimHniniroum filtered Image data and 
the C l a filtered image data. 



CNR ConrpartGon of Different Filters 







T-o 


* S 




t»r:*a&E2 

mai*mlriniBr:a*lj6 




^^OtMOCdQO W«n « " 0.72 

i ■ - 



Flgon 5. The mean CNR measurement 
Of each of 53 vessel segments in an MRA 
processed ustng the G\ filter and the 
cross-section fitter as functions of the 
mean CNR measurement of the same ves- 
sel segments in the same MRA processed 
using the roaxtmunvmlnlmuxn Alter. 



a tor in the maximum-minimum and 
crasB-section filters Is sensitive for Une 
structures. Figure 6b shows a hypotheti- 
cal signal intensity profile of a targe vessel 
& cross the cross-section of the vessel. At 
the center of the vessel, the SD operator 
has a value ore (= 2 x 140-137- 137). 
whereas the operator has a value of 20 
|x 2 X 140 - 130 - 1301. This explains 
why the G£ filter has less signal loss In 
large vessels compared with the maxi- 



haw 
* * ♦ 



„ # *J kf u* 



Figure 6. Hypothetical signal intensity 
profiles df a small vessel, (a), and a large 
teasel, (b). along the cross-sections of the 
vessels. In fa), the vessel center has a sig- 
nal intensity of X30; the background has 
a signal intensity of LOO. la (b), the vessel 
center has a signal Intensity or 140; the 
background has a signal intensity of 100. 



lAUin-uitoljuum and cross- section filters, 
as wd found tn our study. 

The G£ filter is expected to have high 
sensitivity tor vessels with diameter 1M- 1 
gtW points or amaJter. If M in the G* fitter 
ts chosen to be larger than the interpola- 
tion factor. L, the C$ filter is expected to 
have a better sensitivity for larger vessels. 
However, choosing M to be larger than L 
might slightly reduce the sensitivity for 
small vessels. Using a large K value also 
increases the computation time consider- 
ably. Therefore. v« recoanroend using M - 
L and JC = L -1 to avoid this signal loss' 
Tor highly curved email vessels and to keep 
the algorithm computationally efficient. 

• CONCLUSION 

Low vessel- background contrast causes 
some of the amah vessels to lose their vis- 
ibility to MtP images. A cylinder model was 
developed to simulate a three-dimensional 
MR angiogram. A vessel enhancement Al- 
ter, G£, using long-range signal correlation 
was presented. The results of applying the 
04 vessel enhancement filter in three-di- 
mensional MRAs demonstrate that this Mi- 
ter suppresses background signal and 
image noise more effectively compared 
with the rrtaxlmum-niintroum alter and 
cross-section filter ie ported in Du et al. 
(12). One study has shown Chat the Gj fil- 
ter increases the vessel-background CNR 
by a factor of 2.52 cornpared with the max- 
1 mu m-nii nfmu m filter and by a factor of 
3.51 compared with the cross-section fil- 
ter- 
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